use TeO, clear

// Standardize outcomes
foreach x in i_ecole harass_school {
	qui sum `x'
	local `x'm=r(mean)
	sum `x' 
	local sd=r(sd)
	gen `x'st=(`x'-``x'm')/`sd'  
}

//Regressions 
estimates clear


reg harass_schoolst i.birthyear i.religion tmuslim interm if bpl=="1101" &female==1 , cl(religion)
eststo m1

reg i_ecolest i.birthyear i.religion tmuslim interm if bpl=="1101" &female==1 , cl(religion)
eststo m2

esttab m* using "TableA2.csv", star(+ 0.1 * 0.05 ** 0.01 *** 0.001) replace ///
		cells(b(fmt(a3) star) se(par)) stats(N r2)  ///
		keep(interm) 
